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We present results for the heavy quark potential computed in SU(3) from magnetic monopoles and from center 
vortices. The monopoles are identified after fixing SU(3) lattice configurations to the majdmal abelian gauge. The 
center vortices are identified after using an indirect center gauge fixing scheme which we describe for SU(3). Z(3) 
center vortices are extracted and used to compute the potential. The values of the string tensions from monopoles 
and vortices are compared to the full SU(3) string tension. 



A clear understanding of the confinement 
mechanism in QCD is still elusive. Two particu- 
lar mechanisms have gained the most attention: 
magnetic monopoles and center vortices. Here we 
present results of calculations in SU(3) pure gauge 
theory of the heavy quark potential using meth- 
ods which test these two different mechanisms. 



1. Monopoles and the Maximal Abelian 
Gauge 

The maximal abelian gauge is obtained in 
SU{3) gauge theory by maximizing the lattice 
functional 

S^Yl \[Tr{U,{x)XsUl{x)Xs) 

X.fl 

+ Tr{U^{x)XsUlix)Xs)], 

where As and A3 are the diagonal generators of the 
fundamental representation of SU(3). The maxi- 
mization is done by performing local SU{3) gauge 
transformations on the links, which in practice is 
done using an SU{2) subgroup update |l[]. 

Once the maximal abelian gauge is obtained, 
two independent abelian gauge fields are pro- 
jected out, leaving [/(I) x C/(l) gauge fields. Mag- 
netic monopoles are then identified in the abelian 
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configurations. The lattice monopoles are used to 
directly compute the heavy quark potential[^-Q. 



2. Indirect Center Gauge 

The indirect center gauge was used in the 
first lattice studies of projected center vortices 
in S't/(2)||]. The indirect center gauge is ob- 
tained in the following way. A non-abelian lat- 
tice gauge configuration is fixed to the maxi- 
mal abelian gauge. Abelian projection is per- 
formed obtaining the abelian subgroup configu- 
rations: In SU{2) this procedure obtains U{1) 
gauge fields, while in SU{N) it yields [C/(l)]^-^ 
gauge fields. The remaining U{1) configurations 
are then gauge fixed so that the gauge fields are 
brought as close as possible to the elements of the 
center group of the original SU{N) theory. For 
the case of SU{2) this corresponds to maximizing 

i?2 = ^cos^ {0{x,^i)) , 

where 9{x,^) parameterizes the abelian links. 
The maximization is obtained by iteratively per- 
forming U{1) gauge transformations. The final 
step to extract the vortices involves another pro- 
jection of the lattice gauge fields onto the group 
elements of Z{2) £ {— 1,+1} x 1, where 1 is the 
unit matrix. 
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In SU (3) after maximal abelian gauge fixing 
and projection is done, the links are elements of 
t/(l) X U{1). The U{1) X U(l) links can be pa- 
rameterized by the matrix: 

U^{x) = 

U{1) X U{1) gauge transformations are then ap- 
plied to these links to bring them as close as pos- 
sible to the elements of the SU{3) center group, 

Z{3) e {exp(-i27r/3),l,exp(+i27r/3)} x 1. 

The necessary gauge transformations can be 
parameterized using the fundamental, diagonal 
generators of SU{3): 

g{x) = e'^8a2(x)g*A3ai(:c)^ 

where some numerical factors are absorbed into 
the a's and we use 

As = diag(l,l,-2) 

A3 = diag(l,-l,0). 

The center gauge in SU (N) is obtained by max- 
imizing the quantity 

i?jv -^|Tr[/^(x)|'. 

For the indirect case in SU{3), this is ac- 
complished by applying the U{1) x U{1) gauge 
transformations described above to the diagonal 
[/(I) X U{1) links to maximize the quantity 

R3=J2 3 -f 2008(6*1 (x, Ai) -6l2(x,Ai)) 

X.fl 

Gauge transformations must be applied itera- 
tively to each site independently to maximize the 
global quantity. There are two local gauge condi- 
tions that must be met when the center gauge has 



been obtained. The first condition comes from 
the As contribution: 

E^a;[sin(20i(a;,Ai)+02(x,/i)) 

+ sm{2e2ix,fi)+ei{x,fi))] =0, 

where d~ is the lattice derivative 

a;/(x) = /(x)-/(x-A). 

The second gauge condition to be met comes from 
the A3 contribution: 

J:^d-[2sinie,{x,fi)~e2{x,fi)) 
+ sm{29i{x,fi) +e2{x,^)) 

-sm{2e2{x,^l) + ei{x,^))] =o. 

Numerically, neither term can actually be zero. In 
practice, numerical iterations are continued until 
both terms are on the order of 10^^. In the con- 
tinuum these two gauge conditions are equivalent 
to the Landau gauge: 

d^Oi{x,fi) = 
<9p 6*2(2;,^) = 

3. Calculations 

We generated 400 equilibrium SU{3) lattice 
configurations on a 10^ X 16 lattice with /3 = 5.9. 
At this value of (3 the correlation length is about 4 
lattice spacings. Using these configurations and 
smearing techniques]^ we computed the heavy 
quark potential and measured the full SU (3) 
string tension, aj. We found a/ — .068(3). 

Next, the 400 configurations were fixed to the 
maximal abelian gauge and the two species of 
abelian fieldsQ were projected out. Overrelax- 
ation was used with w = l.?]!]. In each abelian 
configuration Dirac strings were located using the 
1x1 plaquettes and the two species of monopole 
configurations were obtained. The monopole con- 
figurations were used to compute the heavy quark 
potential (see Fig. (2) below). Results for both 
monopole species were found to be the same. The 
string tension from the potential computed using 
monopoles was found to be (Jmon — .047(2). This 
is about 30% lower than the full SU{3) string 
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tension tr/. We believe this is the first quantita- 
tive calculation of the SU(3) string tension from 
monopoles, although hints that the monopole 
string tension is low have appeared before [Q. 

After gathering the abelian configurations, we 
next applied our abelian center gauge fixing algo- 
rithm. No closed form solution of the parameters 
ai and a2 could be found for the transcendental 
equations resulting from maximizing R3. The a's 
were determined iteratively. The value of the nor- 
malized abelian trace computed from using the 
C/(l) X U{1) links was observed to change from 
an average value of 0.525 before center gauge fix- 
ing to 0.922 after gauge fixing. 
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Z{3) configurations. The Wilson loop data is pre- 
sented in Figure (1). Fits to the potential com- 
puted using the Z{3) Wilson loops gave a string 
tension of (Tvort = .060(4). Thus, we find the 
string tension from indirect center projection to 
be too low. 
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Figure 2. Comparison of the three heavy quark 
potentials, computed directly from the S'f7(3) 
gauge fields, using monopoles, and using pro- 
jected center vortices. 



Figure 1. Wilson loops computed using Z(3) links 
obtained after indirect center projection. The 
loop values are for i? = 1 to i? = 5 ordered from 
bottom to top, respectively 



From the center gauge fixed configurations, the 
links were replaced by the ^(3) phase each was 
closest to. Wilson loops were computed using the 



In Figure (2) we compare the heavy quark 
potentials. The potentials computed from the 
monopoles and projected center vortices show lit- 
tle Coulomb behavior at small R as compared to 
the full SU (3) potential. This has been observed 
in other calculations in SU(2) and 
Calculations which demonstrate the confinement 
mechanism are not required to reproduce the 
short distance behavior of the potential, only the 
string tension. 
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4. Conclusions 

We have presented results for the string ten- 
sion computed using monopoles and center vor- 
tices. The calculations use operations defined 
at the level of a single lattice spacing. Neither 
monopoles nor center vortices reproduce the full 
SU{3) string tension. There is still a need to re- 
fine the methods used here if one is to discard or 
accept monopoles or center vortices as the con- 
finement mechanism in QCD. We conclude by 
discussing possible reasons to suspect systematic 
errors exist in these calculations. 

Monopoles and vortices are computed using ob- 
jects defined on the single link level. The results 
for the string tension may be corrupted by un- 
physical degrees of freedom due to misidentifica- 
tion of monopoles and center vortices by using 
operators defined at the single link level. Also, we 
may not be using volumes large enough to sup- 
press finite volume effects since monopoles and 
center vortices are objects that likely extend over 
many lattice spacings. To investigate the first 
problem, we need to construct other methods for 
identifying monopoles and center vortices. For 
the second problem we need to work on larger 
lattices. Using larger lattices is difficult due to 
the computationally intensive nature of the gauge 
fixing methods. 

The numerical procedures for the maximal 
abelian and center gauge fixing suffer from Gribov 
copy problems. In particular, for the direct and 
indirect center gauge fixing procedures in SU (2) 
the effects of Gribov copies on the numerical value 
of the string tension are large pT|. 

The effect of taking account of the Gribov am- 
biguity for SU{2) for all cases investigated so 
far, is to decrease the corresponding string ten- 
sion. So for monopoles in the maximal abelian 
gauge, the monopole string tension for one gauge 
fix/configuration is quite close to the full SU{2) 
value H], and decreases by 0(10%) when account 
is taken of the Gribov ambiguity by generating 
many copies/configuration Our result for 

the 5C/(3) monopole string tension is consider- 
ably lower than the full SU (3) value even at one 
gauge fix/configuration, and would presumably 
be reduced somewhat more when gauge ambigui- 



ties are accounted for. Thus the SU{2) and SU (3) 
cases appear to be rather different. Given the 
large noise level of our data, we did not succeed 
in extracting a string tension at the U{1) x U{1) 
level, before projecting to monopoles or vortices. 
Given our low results for the monopole and center 
vortex string tensions, this is now an important 
quantity to measure. A failure of abelian projec- 
tion at the U{1) X U{1) level for SU{3) would have 
important implications for ideas on confinement 
in SU{3). We plan to return to this question in 
future work. 
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